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■^ Abstract 

o 

^-H In a variety of well motivated models, such as two Higgs Doublet Models (2HDMs) 

ff^ and the Minimal Supersymmetric Standard Model (MSSM), there are neutral Higgs 

^^ bosons that have significantly enhanced couplings to b-quarks and tau leptons in com- 

I parison to those of the SM Higgs. These so called non-standard Higgs bosons could be 

1" copiously produced at the LHC in association with b quarks, and subsequently decay 

• ^H into b-quark pairs. However, this production channel suff'ers from large irreducible QCD 

rN backgrounds. We propose a new search strategy for non-standard neutral Higgs bosons 

C^ at the 7 TeV LHC in the 3b's final state topology. We perform a simulation of the signal 

and backgrounds, using state of the art tools and methods for different sets of selection 

cuts, and conclude that neutral Higgs bosons with couplings to b-quarks of about 0.3 or 

larger, and masses up to 400 GeV, could be seen with a luminosity of 30 fb~ . In the 

case of the MSSM we also discuss the complementarity between the 3b channel and the 

inclusive tau pair channel in exploring the supersymmetric parameter space. 



1 Introduction 

The origin of electroweak symmetry breaking, leading to the generation of mass for the quarks, 
leptons and the weak gauge bosons is one of the most outstanding questions in high energy 
physics. In the Standard Model (SM), the spontaneous breakdown of the electroweak sym- 
metry is induced by the introduction of a scalar doublet field that acquires a non-vanishing 
vacuum expectation value. A physical particle appears in association with the Higgs mecha- 
nism, namely the Higgs boson. The tree-level couplings of quarks and leptons to the Higgs 
boson are then well defined and are proportional to the quark and lepton masses and inversely 
proportional to the vacuum expectation value of the Higgs field. Hence, apart from the top 
quark, all quarks and leptons have small couplings to the SM Higgs boson. Therefore the 
production cross section of a SM Higgs in association with all SM fermions, apart from the 
top quark, is too small to be detectable at hadron colliders. 

It is very likely, however, that the electroweak symmetry breaking sector is more compli- 
cated than just a single Higgs doublet. One of the simplest extensions of the SM is the two 
Higgs doublet model (2HDM). In this case, quarks and leptons receive contributions to their 
masses coming from both Higgs doublets. If the vacuum expectation value of one of the Higgs 
doublets is small, its coupling to some of the quarks can be very large. Flavor physics puts 
however additional constraints on these extended Higgs sectors: in order to suppress large 
flavor- changing neutral-current (FCNC) interactions, either the coupling of one of the two 
Higgs doublets to fermions with a given electric charge is suppressed, or there is an alignment 
between the couplings of the fermions to the two Higgs doublets. The so-called Type II 2HDM 
belongs to the first class of models: up-type quarks and neutrinos couple to one Higgs doublet, 
iJ„, and down-type quarks and charged leptons couple to the other, H^. 

The minimal supersymmetric extension of the SM (MSSM) contains two Higgs doublets. 
At the tree-level, the MSSM Higgs sector is a Type II 2HDM. However, fermion couplings to 
both Higgs doublets are induced at the loop level. Once Supersymmetry (SUSY) is broken, 
dangerous FCNC interactions are generated, but since they are proportional to loop-induced 
couplings, they tend to be suppressed. Due to the large top quark mass, the vacuum expecta- 
tion value of the Higgs field that only couples to the up sector at the tree level {Hy) cannot be 
much smaller than the SM one. Defining tan/3 as the ratio of the vacuum expectation values 
of the two Higgs fields (tan/3 = {Hu)/{Hd)), this implies that tan/3 should be of the order of 
or larger than one. 

In Type II 2HDMs as well as in supersymmetric models, large values of tan /3 imply a large 
coupling of the b quark to the non-standard Higgs bosons, resulting in both a large production 
cross section of Higgs bosons in association with b quarks as well as a large branching fraction 
of the Higgs bosons decaying into b quarks. In the supersymmetric case, the precise value of 
the coupling depends not only on tan/3, but also on SUSY-breaking effects. These effects can 
modify both the b and r couplings to the non-standard Higgs bosons, and hence the branching 
ratio of these Higgs bosons decaying into b quarks and r leptons. 

Non-standard Higgs boson production at the LHC has been mainly studied through inclu- 
sive Higgs boson decays into r leptons, since this channel has a reasonable signal-to-background 
ratio. Currently, the LHC experiments are setting strong bounds [1, 2] on light non-standard 
neutral Higgs bosons at moderate or large values of tan /3, surpassing the previous bounds [3, 4] 
set by the Tevatron experiments. In this work, we shall study the associated production of 



a non-standard neutral Higgs boson with b quarks at the LHC, with the Higgs boson subse- 
quently decaying into b quarks. The process involves the production of at least three b quarks in 
the final state and, for large values of tan /3, the production cross section may be sizable. This 
search channel suffers from a large irreducible background, since the QCD bb+X production 
cross section is much larger than the one associated with Higgs production. Hence, previ- 
ous experimental studies of this channel have been mainly performed at the Tevatron [5, 6], 
where the backgrounds are easier to control, but no CMS or ATLAS analysis is available at 
present. This search channel is challenging because the b quark produced in association with 
the Higgs boson typically has low transverse momentum {pt) and triggering on such soft b 
jets, especially for low Higgs boson masses, is difficult at high instantaneous luminosity. After 
satisfying the trigger requirements, demanding that two of the b jets reconstruct the Higgs bo- 
son invariant mass helps to improve the signal significance. However systematic uncertainties 
can still be an issue due to the small signal-to-background ratio. Previous theoretical studies 
focused on the prospects for the discovery of a non-standard neutral Higgs boson in the 36 
and 46 channels at the 14 TeV LHC [7, 8, 9]. In particular supersymmetric models such as 
anomaly-mediated supersymmetry breaking (AMSB) models and gauge-mediated supersym- 
metry breaking (GMSB) models have been studied [9]. The aim of our paper is to analyze 
the reach at the 7 TeV LHC for Higgs bosons arising in generic 2HDMs. In particular, we 
determine the required effective coupling of the neutral Higgs boson to b quarks to have a 
possible discovery at the 7 TeV LHC with 30 fb^^ of data. We also consider the specific case 
of the MSSM and investigate the complementarity to the A ^ rf searches. 

In section 2, we shall present the necessary theoretical background and emphasize the 
differences between the several two Higgs doublet extensions discussed above. In section 3 
we shall study the reach for non-standard neutral Higgs bosons at the 7 TeV LHC in the 36 
channel. We describe our simulation of signal and background and the proposed selection 
cuts, and then discuss the expected reach in specific 2HDMs, as well as in the MSSM. We 
reserve section 4 for our conclusions and outlook. 

2 The 36 channel in 2HDMs 

In a 2HDM the most generic Yukawa couplings of the two Higgs doublets with SM quarks and 
leptons can be written as 

-CYuk = Vu H^QU + Vd HdQD + y„ hIQU + m HlQD + y, H^LE + y, HILE + h.c. , (1) 

in which H^ and Hd are the two Higgs doublets with hypercharge 1/2 and -1/2, respectively. 

A generic structure of the four Yukawa couplings leads to Higgs-mediated FCNC interac- 
tions already at the tree level. However New Physics (NP) effects in fiavor transitions can be 
reduced by imposing the alignment of up, down and lepton Yukawa couplings [10] or, more 
generically, the Minimal Flavor Violation principle [11]. 

We shall introduce the variables e/, parameterizing the relation between the yf and yj 
couplings in alignment models, 

yt = ttyu Vb = tbyd, tJT = tryr , (2) 

with generic fiavor independent et^b,T coefficients. In the following we will focus on the couplings 
of the three neutral Higgs bosons with the third-generation down-type quarks and leptons. 



Assuming that there is no CP violation in the Higgs sector, the couphngs of the neutral Higgs 
bosons with b quarks are given by [12] 
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where rhb is the running b quark mass and a is the mixing angle between the two scalars h, H. 
The corresponding couplings of the Higgs bosons with the third-generation charged leptons 
are obtained with the simple exchange b ^ t. Hence, in generic aligned 2HDMs the couplings 
of the pseudoscalar Higgs boson with b quarks and r leptons can be parametrized by two 
independent effective couplings tan (3^^ and tan (3^^ 

9 .^„ tan/3 f _^ e^ \ , , ,j ,, _ g ^ _„ ^b 



g _ tan/3 / ^t \ , . ^ •, 9 - 

2MH;^^l + e.tan/3 [' ' t^J ^^^^^^^^ " ^^-.tan/3,,(Azf75r) . (5) 

Moreover, in the decoupling limit, arising at large values of rriA and tan/3, cos a ~ sin/3, 
sin a ~ — cos /3 and consequently the coupling of the heavy scalar H with b quarks (r leptons) 
is also governed by tan/3^g (tan/3jg). The coupling of the light scalar h is instead SM-like in 
this limit. 

The total production rate of b quarks and r pairs mediated by the production of a CP-odd 
Higgs boson (as well as by the heaviest CP-even Higgs scalar) in the large tan/3 regime can 
be approximated by [13] 

ftan /3^ Vfh'^N 
a{bb ^ A)BniA -^ bb) ~ a{bbh)sM{i^nPl^) \^ oX - u \b\.-2.r ' (6) 



a{gg,bb ^ A)Bn{A ^ tt) ~ a{gg,bb ^ /i)sM(tan/3,^ff) %^ _^ ^^ ^^^^^ '^;^^_^ ^ (7) 



(tan/3^-g)2m2 + (tan/3^^g)2m2iV, 

,b ^2 (tan/3;ff)2m^ 

(tan/3-jj)2m2 + (tan/3,^g)2m2Ar, 



where Nc is the number of colors {Nc = 3) and <j{bbh)sM and a{gg,bb — )■ /i)sm denote the 
values of the corresponding SM Higgs boson production cross sections for a Higgs boson of 
equal mass^. 

The MSSM at the tree level is a particular 2HDM of Type II, however at the one-loop 
level also the Yukawa couplings yu,yd,ye are generated. In particular the dominant threshold 
corrections to the b quark mass are arising from gluino-sbottom one-loop diagrams and from 



^In Eq. (7) we have neglected the contribution to the production cross section coming from the top-quark 
loop diagram. The corrections arising from the interference terms between the top-quark and 6-quark loop 
diagrams amount only to a few percent [13]. 



chargino-stop loops, resulting in e^ ~ eo + y'^ey with [14]-[18] 
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in which Mg is the gluino mass, M^, and M^. the sbottom and stop masses and /i the Higgsino 
mass parameter. At and yt are the top trilinear term and the top Yukawa coupling, respectively. 
Similarly, the corrections to the r mass are dominated by wino and bino exchange con- 
tributions that are usually small, since they are suppressed by the electroweak coupling and 
have the form [19, 20] 

e^~^^M2/(M,^,M2,/i), (10) 

with M2 the wino mass, M^^ the sneutrino mass. 

The effective couplings of the CP-odd and heavier CP-even Higgs bosons of the MSSM 
with h quarks and r leptons are then given by Eqs. (4) and (5) with the resummation factors 
£{, and e^ given just above. As a result, contrary to generic aligned 2HDMs, in the MSSM the 
two couplings will be uniquely determined, once the supersymmetric spectrum is specified. 

3 Early LHC prospects for the 36 channel 

3.1 Simulation of Signal and Background 

Signal and background processes are modeled using the MadEventS [21] event generator in- 
terfaced with Pythia 6.4 [22] for parton showering and hadronization, using a matrix-element 
parton-shower matching algorithm to avoid double-counting of partonic configurations. We 
allowed for up to two additional partons in the final state of the hard process when following 
the shower-ZcT scheme outlined in Refs. [23, 24], with /c^- matching scale of 30 GeV. These sam- 
ples were generated for pp collisions at ^/s = 7 TeV using the CTEQ6L1 parton distribution 
functions (PDFs) [25]. 

The QCD production of multiple heavy quarks is the main source of background. We 
generated two separate QCD background samples: a bbj + X {j = u, d, s, c, g) sample in which 
the additional partons used in matching to the parton-shower are light or charm quarks and 
a "36" {bbb + X and bbb + X) sample in which the additional partons could be either light or 
heavy quarks.^ 

We cluster particle jets using the anti-Zcr algorithm implemented in FastJet-2.4.3 [26] with 
a radius parameter R = 0.4. To better simulate the experimental b jet energy resolution, we ex- 
clude neutrinos from jet clustering, which in the case of semileptonic b decays can carry away a 
significant fraction of energy. Furthermore, we apply a jet energy smearing of 100% / ^^/E/GeV 
to model the typical calorimeter energy resolution of LHC experiments. 



^This separation of QCD background into the bbj and 36 samples does not model b jets with pT below 
^ 40 GeV very well. However once we impose px cuts on the jets as described in the following, the effects are 
only at the ^ 10% level. 
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Figure 1: Invariant mass distributions for each of the 6-jet pairs for events satisfying Selection I and 
assuming an integrated luminosity of 30 fb~^ in pp collisions at -y/s = 7 TeV. The expected distribution 
for the total background (red histogram) is compared to that for a signal with rriA = 200 GeV (blue 
histogram) and niA = 300 GeV (green histogram) assuming tan/3 = 30. The signal expectation has 
been scaled by a factor of ten for visibility. 



Since the signal typically contains three b quarks, flavor tagging becomes an effective 
tool to suppress QCD multijet backgrounds. We assume a constant 6-tagging efficiency of 
60%, a c-jet mis-tag rate of 10% and a light-jet mis-tag rate of 1% [27]. This choice can be 
considered conservative, as LHC experiments have already developed sophisticated 6-tagging 
algorithms [28] exceeding the performance assumed in this paper. The low mis-tag rate of c- 
and light-jets leads to the bbj and 36 backgrounds being comparable once three 6-tagged jets 
are required. 

We consider two sets of event selection criteria: 

1. Selection I: events are required to have exactly three 6-tagged jets with px > 60 GeV 
and \ri\ < 2.0. 

2. Selection II: events are required to have exactly three 6-tagged jets with px > 50 GeV 
and \ri\ < 2.0, and the leading 6-tagged jet to have pt > 130 GeV. 

In both cases jets are required to be relative central to ensure they are contained within 
the tracker volume and can therefore be tagged with high efficiency [27]. 

The high instantaneous luminosities delivered by the LHC has forced to raise threshold in 
the trigger menus in ATLAS and CMS. As a result, the ffist selection is somewhat optimistic 
in that the jet pt requirements may be too low for these events to satisfy trigger requirements 
with high efficiency. For instance, events satisfying Selection I would have a low efficiency to 
satisfy the 26/2j ATLAS trigger [29] requirements, which appears as one of the most suitable 
unprescaled triggers for this topology. On the other hand. Selection II would be more rep- 
resentative of the kind of minimum jet pt requirements applied by the 26/2J ATLAS trigger 
once trigger turn-on effects are considered. Nevertheless, we consider Selection I to explore 
the potential sensitivity gains at low ttia values, which could motivate designing an optimized 
trigger strategy for such lower -px events at the LHC experiments. 

In most part of the parameter space under consideration, there is no obvious way to 
accurately distinguish the pair of b jets coming from the decay of the heavy scalar and the b 
jet produced in association with it. However Selection II with masses mA ^ 260 GeV represents 
an exception: the b jet produced in association with the Higgs boson is the leading one (6i) 
because the b jets coming from the decay of the Higgs are too soft to satisfy the selection 



criteria. Therefore in Selection II for low rriA the invariant mass of the second-leading b jet (62) 
and the third-leading 6 jet (63) reconstructs the Higgs mass, while for large rriA (^a ~ 260 GeV) 
the invariant mass of 61 with 62 or 63 reconstructs the Higgs mass. Nevertheless to improve 
the acceptance, we consider all three possible combinations of 6-jet pairs and require that 
the invariant mass of at least one of them is within a window around the peak of the Higgs 
boson invariant mass distribution. Figure 1 compares the invariant mass distribution between 
signal and background for each of the possible 6-jet pairs. As it can be appreciated, the signal 
distribution consistently peaks at values lower than the physical mass of the Higgs boson due to 
a combination of the jet smearing^ and energy loss via neutrinos, with the PDF suppression 
involved in producing such a massive resonances. This effect becomes more significant for 
heavier resonance. If an excess in the invariant mass distribution of a pair of b quarks is 
discovered at the LHC, one needs to extract the actual mass of the resonance through proper 
simulation. The chosen central value of the mass window is shown in Table 1 for each of the 
niA values considered. 

Additionally, we studied the effect on varying the width of the mass window about the 
peak and found that a typical width of \Mbb — "nipeakl < 25 GeV yields good results across the 
whole mass range. Increasing the width of the mass window for heavier masses and reducing 
the width for lower masses can lead to a improvement in the significance of only of a few 
percent. Finally, we checked that imposing a cut on AR between any two b jets did not lead 
to a marked improvement in the signal significance. 



niA (GeV) 


150 200 250 300 400 500 


"^peak (GeV) 


150 190 230 250 350 450 



Table 1: Physical masses, itla, used in our signal samples, and the corresponding central 
values of the mass window. The mass window used in the signal selection is chosen to be 
\Mbb - mpeaki < 25 GeV. 



3.2 Prospects and Significance 

As discussed in Sec. 2, in a generic 2HDM the couplings of the A, H bosons with b quarks 
and with r leptons are independent and parametrized by the effective couplings tan/3^g and 
tan/3Jfj, respectively. As shown by Eq. (6), the cross section times branching ratio for non- 
standard neutral Higgs bosons produced in association with a b quark which subsequently 
decay into a pair of b quarks has only a mild dependence on the choice of tan/J^Tfj. Without 
lost of generality we fix tan P^^ = 5 to be in agreement with the present bounds coming from 
LHC A,H ^ Tf searches [1, 2]. 

In addition, we are focusing on the parameter region with sizable tan/J^g and moderate 
values of rriA- In this case, the heavy CP-even scalar H and the pseudoscalar A can only 
be slightly split in mass, so that the two Higgs bosons will appear at the LHC in the same 
resonance region with combined cross section. Hence, for our analysis, the only two relevant 

■^We checked that, in the range of tan /3 we are considering for our analysis, the effects on the b-pair invariant 
mass distributions coming from the physical width of the Higgs are negligible, if compared to the width of the 
invariant mass distribution coming from jet smearing. 



free parameters are rriA ~ ttlh and tan/J^g. To compute the rate of the signal, we double the 
cross section obtained for the CP-odd Higgs and use the narrow width approximation which 
is valid in the entire mass range we consider as long as tan/J^g is not too large (tan/3^g < 80). 
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Table 2: Expected number of background and signal (at tan/3^g = 30) events per 30 fb^ of data at 
the 7 TeV LHC, after imposing Selection I presented in the text (above double line) and after the 
mass window selection presented in Table 1 (below double line). The first row shows the total events 
in each channel before event selection criteria are imposed. 





36 


bbj 




S 


gnal (mA in 


GcV) 










150 


200 


250 


300 


400 500 


After matching /lO"^ 


4800 


2.2 X 10^ 


420 


180 


90 


45 


14 5 


Selection II 


24000 


42000 


1200 


1650 


2100 


1850 


850 120 


m^ = 150 GeV 


6300 


11000 


800 










niA ^ 200 GeV 


10000 


19000 




1350 








niA = 250 GeV 


12000 


20500 






1700 






niA ^ 300 GeV 


11000 


20000 








1200 




niA = 400 GeV 


4800 


9000 










390 


TTiA ^ 500 GeV 


1900 


2900 










45 



Table 3: Expected number of background and signal (at tan/3gg = 30) events per 30 fb^ of data at 
the 7 TeV LHC, after imposing Selection II presented in the text (above double line) and after the 
mass window selection presented in Table 1 (below double line). The first row shows the total events 
in each channel before event selection criteria are imposed. 

In Tables 2 and 3 we present the number of signal and background events per 30 fb~^ in 
each of the test mass windows, assuming tan/^g^ = 30. In Table 4, we compare the signal 
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Table 4: Signal (at tan /3^jj = 30) to background ratio and significance S/\/B per 30 fb of data at 
the 7 TeV LHC, using the two Selections presented in the text. 



statistical local significance^ from applying Selection 1 with that of Selection 11, assuming an 
integrated luminosity of 30 fb~^. We can see that, as expected. Selection 1 has a markedly 
better statistical sensitivity for m^ < 300 GeV. 



Using Eq. (6), it is straightforward to generalize these results to different values of tan /3, 



eff- 



In Fig. 2 we present the accepted signal cross section, after that all cuts are implemented, 
and the statistical significance at 30 fb~^ LHC as a function of the two free parameters of the 
theory, tan (31^ and m^ for both Selection 1 and Selection 11. The results are encouraging. For 
a 7 TeV LHC run with a total integrated luminosity of 30 fb~^, we can probe a large parameter 
region of the 2HDM by searching for heavy Higgs scalars in the 36 final state. For example, 
applying Selection 1, a coupling of the pseudoscalar with h quarks of the order ~ 0.3, and 
hence tan/J^g ~ 30, could lead to a 5(T significance for a Higgs boson mass up to ~ 370 GeV, 
with 30 fb^^ of data. However, from Table 4, we also see that for moderate tan/3gg, the final 
S/B < 0.1. Therefore the potential systematic uncertainties which have not been accounted 
for in our analysis could make signal identifications challenging. We expect that a detailed 
experimental analysis exploiting the sideband regions in the invariant mass distribution to 
constrain systematic uncertainties in the background prediction can nevertheless achieve a 
high sensitivity. In any case, for larger values of tan/3^g(~ 60), and hence a coupling ~ 0.6, 
we can achieve a significance ^ 10 with 30 fb^^ of data in almost the entire range of masses 
considered. For such large effective couplings also the signal-to-background ratio would be 
more favorable, S/B < 1/5. 

Comparing the first and second rows of Fig. 2, we note that there is a difference in shape of 
the contours. Applying Selection 11, the best reach is for tua ~ 250 GeV. On the other hand, 
the reach when using Selection 1 monotonically decreases as the value of tua increases, which 
is the result of the rapid decrease in signal rate. This difference in the shapes of the accepted 
signal cross section and of the statistical significances are related to the different cuts on the 
three highest pt jets. First of all, the h jet produced in association with the Higgs boson has 
a steeply falling distribution, suppressed by the PDFs. Therefore, it is unlikely that this jet 



^Note that in Tabic 4 and in the rest of the paper we arc only presenting the local signihcance for a Higgs 
with mass in one of the selected mass windows. The study of the lookclscwherc effect goes beyond the scope 
of this paper. 
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Figure 2: Accepted signal cross section in pb (left column) and statistical significance S/vB (right 
column) in a generic 2HDM (assuming tan (3'^^ = 5 ) at the 7 TeV LHC for an integrated luminosity 
of 30 fb~ . In the top row, we show the results after imposing Selection I and the mass window cuts 
detailed in Table 1. In the bottom row, the results from applying Selection II and the mass window 
cuts presented in Table 1. 



can satisfy the cut of px > 130 GeV on the leading jet in Selection II. This explains why 
Selection I leads to better sensitivity, in particular in the low mass region. On the other hand, 
for iJiA ~ 250 GeV, it is easier for the b jets from Higgs decay to be the leading jet and to 
satisfy this cut. This explains why the reach in Selection II is better for higher Higgs masses 
than for lower ones. This effect suggests that an asymmetric jet energy cut, similar to that of 



Selection II, could have advantages. The pt cut on the leading jet could be optimized further, 
such as requiring it to be proportional to the target signal mass. 

The possibility of detecting the pseudoscalar and the heavy scalar of the MSSM in the 36 
channel deserves a special discussion, since in the MSSM the effective couplings tan/J^g and 
tan/3jg defined in Eqs. (4) and (5) are determined, once the SUSY spectrum is specified. 

As discussed in Sec. 2, contrary to generic 2HDMs, in the MSSM the coupling of the 
pseudoscalar Higgs with b quarks and r leptons depend equally on tan /3 but have a different 
dependence on corrections arising at the one-loop level. Typically, for gluinos at the TeV 
scale, stops, sbottoms and charginos at a few hundred GeV and At of the order 1-2 TeV, e^ is 
at the few % level. On the other hand, in the lepton sector typically e,- ~ (9(10~^).^ 

In our numerical analysis, we choose two representative scenarios: the first with e^ = 
and eb = —1/60, and the second with e^ = and ej, = —1/30. Both scenarios can be achieved 
in models with a large and negative fi term (see Eqs. (8) and (9)). The effects of introducing 
a small but non-zero e^ will not significantly modify our conclusions. These scenarios are 
presented in Fig. 3. The plots on the left represent the case e,- = and ej, = —1/60; the ones 
on the right e,- = and eb = —1/30. In white we present the bound on tan/3 coming from the 
requirement that the narrow width approximation is valid [Ta ^ i^)- 

It is interesting to compare the shape of the exclusion bound from the LHC A ^ rf search 
to the one of the constant significance contours for the A ^ bb channel. For large values of rriA, 
the present CMS bound [2] obtained with ~ 5 fb^^ of data is weaker than what was expected. 
As a consequence, the slope of the present rf exclusion bound (solid red line in the figure) is 
much steeper than the slope of the A ^ bb constant significance contours, at large values of 
rriA- Differently, the slope of the expected bound projected at 30 fb~^ (dashed red line in the 
figure) gets much closer to the slope of the A ^ bb constant significance contours, especially 
at small values of tb (e^ = —1/60). For larger e^ (e^ = —1/30) one can still note a difference in 
the slopes, that is mainly due to the different scaling in tan /3 of pp — )■ bbA with A ^ bb and 
pp — )■ bbA with A — )■ rf, as shown by Eqs. (6), (7) once tan/3^jj and tan (3^^ are replaced with 
their expressions in (4), (5) and e^ is set to 0. Thanks to this different scaling, the 3b channel 
can be used to probe heavier pseudoscalar masses than the rf channel. Whether the 2a excess 
at high mass from CMS turns out to be a hint for a signal or just the result of a statistical 
fluctuation, it seems imperative to probe this mass range with a channel with comparable or 
better sensitivity, able to provide complementary information on the MSSM preferred region 
of parameter space. 

4 Conclusions 

In this article we have studied the associated production of non-standard neutral Higgs bosons 
with b quarks at the LHC. Considering the Higgs boson decay into b quarks, we have analyzed 
the reach of the 7 TeV LHC collider. After applying a rather simple event selection cuts, a 
manageable signal-to-background ratio could be obtained, helped by the enhanced production 
cross section obtained at large values of tan /3, which allows for a statistically meaningful reach 



^Notice, however, that scenarios with Kght third-generation sleptons and large values of the /i parameter 
can also lead to values of e,- at the few % level [30] . 



10 



ej=-l/60, Selection I 



40 ri 



35 
30 

25 

5 20 

15 

10 






ej=-l/30, Selection I 


?s 






^^^^^^HH^^^ 


20 

QQ. 15 


% 


10 




10 


"^1^^,,^^^^ 


■ 


5 


fc 


^j 




, 1 



150 200 250 300 350 400 450 500 
rtiA [GeV] 



150 200 250 300 350 400 450 500 
niA [GeV] 



e6=-l/60. Selection H 



e6=-l/30. Selection H 



=§• 



40 r; 



35 
30 

25 
20 
15 
10 




150 200 250 300 350 400 450 500 
niA [GeV] 



25 



20 



QO 15 



10 




150 200 250 300 350 400 450 500 
niA [GeV] 



Figure 3: Statistical significance at the 7 TeV LHC for an integrated luminosity of 30 fb~^ in two 
different scenarios: e^ = —1/60 (left panels) and e;, = —1/30 (right panels). The results for Selection 
I and II are shown in the top and bottom rows, respectively. The red solid (dashed) line represents 
the present (projected at 30 fb~^) bound on non-standard Higgs bosons decaying to rf [2]. The area 
below the white solid line corresponds to the region of validity of the narrow-width approximation 
(r^ < m^/lO). 

at moderate values of the LHC luminosity. In particular, we find that Higgs bosons with a 
coupling to h quarks of about 0.3 or larger (namely tan/3gg = 30 or larger), and with a mass 
up to 400 GeV could be discovered with a luminosity of 30 fb~^. We expect that the run at 
8 TeV will enhance the reach by about 10-15%, although a precise estimation of the reach 
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depends on many details which are beyond the scope of this study. 

We have studied the discovery potential using two different sets of cuts. In general, the b jet 
produced in association with the Higgs boson tends to be soft, driven by the suppression from 
the steeply falling PDFs. At the same time, the pt of the b jets from the Higgs decay is closely 
correlated to the mass of the Higgs boson. Therefore, in particular in the low mass region, a 
somewhat lower threshold on the total jet pt will enhance the discovery reach. At the same 
time, it could be beneficial to use an asymmetric pt selection criteria with the requirement 
that the cut on hardest 6-jet px is correlated with the target Higgs mass. This effect should 
be more prominent for higher Higgs masses. 

We have also studied the discovery potential in SUSY-like scenarios. In this case, the 
corrections to the Yukawa couplings arise at loop level, and there is a correlation between the 
bb and rf search channels. We found that the 36 channel can be important in probing super- 
symmetric scenarios in which SUSY-breaking effects can significantly modify the couplings of 
non-standard neutral Higgs bosons to b quarks and r leptons. In particular we showed that 
the rf channel still has a better reach for lower Higgs boson masses, but the bb channel can be 
used to probe heavier pseudoscalar masses than the tt channel. Furthermore, the 36 channel 
provides an important probe into the coupling of the Higgs boson to b quarks and hence it is 
complementary to the rf channel. 
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